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ABSTRACT 

After reviewing some of the mathematical foundations and numerical difficulties fac- 
ing lattice QCD, I review the status of several calculations relevant to experimental 
high-energy physics. The topics considered are moments of structure functions, 
which may prove relevant to search for new phenomena at the LHC, and several as- 
pects of flavor physics, which are relevant to understanding CP and flavor violation. 
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1 Introduction 



Several areas of research in elementary particle physics (as well as nuclear physics and 
astrophysics) require information about the long-distance nature of quantum chro- 
modynamics (QCD). Sometimes this information can be gleaned from experiment, 
but often what one needs, in practice, are ab initio calculations of the properties of 
hadrons. In some cases one aims for a detailed understanding of QCD in its own 
right. In others, one simply requires a reliable calculation of hadronic properties, so 
that one can study electroweak interactions or new phenomena at short distances. 



including QCD, is to start with a space-time lattice. The spacing between sites 
is usually called a. If the finite grid has Ng x N4 sites, then one has a finite box 
size, L = Nsa, and finite extent in time, L4 = N^a. Quarks are described by lattice 
fermion fields located at the sites, denoted ip(x); gluons are described by lattice gauge 
fields located on the links from x to x + a/2, denoted U^{x). The key advantage to 
the lattice is that local gauge invariance is simple. The fields transform as 



so it is easy to devise gauge invariant actions, i.e., independent of g{x) M. If one 
imagines a smooth underlying gauge potential A^x) (as is used in continuum QCD), 
the relation to the lattice gauge field is 



Continuum QCD is defined from lattice QCD by taking a — > with L and hadron 
masses fixed. Then one takes the infinite volume limit, L — > 00. If one is interested 
in the chiral limit, m q — > 0, it should be taken last. These limits are nothing radical: 
the lattice provides an ultraviolet cutoff, and the finite volume an infrared cutoff. 



The existence of these limits has not been proven rigorously, but, because of 



asymptotic freedom, there is not much doubt that this procedure works. (If not, why 
does QCD work at all?) The lattice formulation makes field theory mathematically 
similar to statistical mechanics and, consequently, provides new tools. For example, 
a finite lattice makes it possible to integrate the functional integral by Monte Carlo 
methods. The expectation value of observable $ may be written 



Mathematical physicists tell us that the best way to define gauge theories, 



ip(x) I— > g(x)ijj(x), i/j(x) I— > ip(x)g (x), 
Ufiix) h-> g(x)U fl (x)g~ 1 (x + afl), 



(1) 
(2) 




(3) 



(4) 
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where Z is chosen so that (1) = 1, and is an abbreviation for all fields, ip, ip, and U. 
The right-most expression is a numerical approximation, with the sum running over 
some set of field configurations. This numerical technique, though only one facet of 
lattice gauge theory, is what most particle physicists mean by "lattice QCD." So, 
this talk is about tools need to make numerical calculations more reliable, and the 
progress being made in calculations needed to interpret "physics in collision." 

It is not so easy to descend from the mathematical high ground down to 
realistic, practical numerical calculations. Difficulties arise because QCD is a multi- 
scale problem. Nature has not only the characteristic scale of QCD, Aqcd, but also 
a wide range of quark masses, leading to a hierarchy 

m q < Aqcd < rn Q . (5) 

As a dynamical scale, rather than a parameter, there is a range for the QCD scale. 
Some benchmarks include the scale in the running coupling A^g ps 250 MeV, typical 
hadron masses like m p = 770 MeV, and the scale of chiral symmetry breaking 
m 2 K /m s = 2500 MeV. The strange quark, with m s ~ 100 MeV, is light, and the up 
and down quarks, with rh = + m u ) = m s /24 and > m u > 0, are especially 
light. The bottom quark, with mj = 4.25 GeV, is heavy, and the top quark, with 
m t = 175 GeV, is especially heavy. One can argue whether the charmed quark, with 
2m c = 2.5 GeV, is heavy or not. 

Cutoffs are needed to put the problem on a computer, and they introduce 
two more scales. The idealized hierarchy is now 

L- 1 «m ? « A QCD < m Q < a -1 . (6) 

It is impractical to expect a huge separation of scales in computational physics. To 
explain why, some simple scaling laws are helpful. The memory required grows like 
NgN 4 = L' i L i /a A 1 and these large exponents come because we live in 3 + 1 space- 
time dimensions. The CPU time needed to update gauge fields in the Monte Carlo 
scales like a~( 4+z \ where z = 1 or 2, and the 4 again comes from the dimension of 
space-time. The CPU time needed to compute quark propagators scales like m~ p 
where p = 1-3. The exponents z and p are non-zero because of properties of our 
numerical algorithms, and it seems difficult to reduce them. 

The first consequence of these scaling laws is that a -1 can be larger, but 
not much larger, than Aqcd- Similarly, m q and L~ 1 can be smaller, but not much 
smaller. A more important consequence is that improved methodology pays off 
enormously. For example, B physics with a" 1 ~ 3Aq C d instead of 3m b saves a 
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factor (mft/AQ CD ) 6 > 2 12 in computing. Such improvements are not attained by 
CPU power, but through new ideas. Moreover, in computational physics the need 
for creativity means that computing facilities must be flexible, not just big. 

So, we see that finite computer resources force us to consider the hierarchy 

L' 1 < m q < A QCD < m Q ~ a" 1 , (7) 

instead of the idealized one. One should emphasize that we know how to get from 
the practical hierarchy (|7|) to physical results. The central idea is to let the computer 
work on dynamics at the scale Aqcd, and to use effective field theories to get the 
rest. The computer runs, by necessity, with finite cutoffs and artificial quark masses. 
With effective field theories, one can strip off the artifice, and replace it with the 
real world. In doing so, one introduces theoretical uncertainties, but effective field 
theories control the error analysis. 

A more thorough exposition of this line of thinking can be found in a recent 
review ||. Here let us emphasize the role of effective field theories by listing some 
of the big ideas in lattice of the last several years: 

• static limit and lattice NRQCD to treat heavy quarks 

• understanding lattice perturbation theory (to match at short distances) 

• non-perturbative implementation of the Symanzik effective field theory 

• continuum HQET to control heavy-quark discretization effects 

• novel applications of chiral perturbation theory 

• understanding chiral symmetry in lattice gauge theory 

All but the last explicitly bring in effective field theories, and it resonates with the 
usage of chiral perturbation theory to extrapolate light quark masses. 

The exception to the rule of effective field theory is something called the 
quenched approximation. Quenched QCD is a model, so the associated uncertainty 
is difficult to estimate. For example, unquenched (usually rif = 2 not 3) calculations 
of hadron masses, decay constants, etc., suggest changes of 0-20%. Fortunately, 
the quenched approximation is going away. Within a few years, I imagine that 
quenched calculations will no longer play an important role in our thinking about 
non-perturbative QCD. 

The rest of this paper looks at some calculations needed to interpret current 
and future experiments. Moments of structure functions can help obtain better 
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parton densities and, hence, better predictions of cross sections at the Tevatron and 
LHC. Lattice calculations of these moments are covered in Sec. [| In flavor physics, 
the focus of many experiments, there is a pressing need for the hadronic matrix 
elements needed in leptonic and semi-leptonic decays, and neutral meson mixing. 
There are many of these, and Sec. |3| covers only a subset: form factors Tb^d* (1) and 
f+(E) to obtain \V c b\ and \V u b\ from B — > D*lv and B — > nlu; and matrix elements 
for neutral B, B s , and K oscillations, need to constrain V t d- Section [| concludes 
with a summary and some thoughts about the future. 

2 Moments of Structure Functions 

The rate for deeply inelastic Ip scattering (DIS) depends on several structure func- 
tions, which we shall generically denote F(x). x is the momentum fraction of the 
struck parton, < x < 1. In perturbative QCD, the Fs can be related to process- 
independent parton densities, which are used to predict cross sections for pp and pp 
collisions. The DIS data peter out for x > x max ~ 0.7, so, as sketched in Fig. |l](a), 
the uncertainty explodes for the highest values of x. High x partons are needed to 
produce high-mass particles. Modern methods for parton densities directly reflect 
the lack of information ||. They need independent (but QCD-based) knowledge of 
the moments to constrain F(x) for x > x max . 

The operator product expansion relates the moments to local operators, 

( l dxx n ~ l F(x,Q 2 ) = C n (Q 2 /fi 2 ) ( P \O n \p)(ri, (8) 
Jo 

where C n is a short- distance Wilson coefficient that is calculated in perturbation 
theory. The matrix element on the right-hand side should be "easy" to calculate: just 
find a lattice operator O nlat with the same quantum numbers as O n and calculate the 




Figure 1: (a) Sketch of the uncertainty in measured structure functions, as a function 
of x. (b) Comparison of chiral extrapolations for (x) u -d; from Ref. j7Z|/ 
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matrix element at several different lattice spacings. We have a detailed description 
of lattice-spacing effects, namely, 



(p\O n \M(a) = Zj n {^a){p\O m \p)(^) + aK nj {p\0' j \p) 

+ aK aF Z~^ J d 4 x (p\TO m qaFq(x)\p) + 0(a 2 ), (9) 

based on an effective field theory introduced by Symanzik two decades ago Q. 
Differences between lattice gauge theory, with a ^ 0, and continuum QCD arise at 
short distances comparable to a. They are lumped into short- distance coefficients 
Z nm , aK n j, and aK a p. The operators on the right-hand side are, in the Symanzik 
effective field theory, defined with a continuum renormalization. Thus, in addition 
to calculating the left-hand side, to obtain (p\O m \p) one must also compute the 
normalization factor Z nm and cope with the terms of order a. 

We know how to calculate Z nm , K a p, and K n j. They depend on details of 
the lattice Lagrangian and the lattice operator O nlat . One can introduce parameters 
(call them csw an d c n j) that directly influence K aF and K n j. Thus, one can ad- 
justed csw and c n j until K a p « and K n j « 0. This procedure is called Symanzik 
improvement |J. There are two ways to compute the Zs and Ks. One is renor- 
malized perturbation theory, which works because the Zs and Ks are short-distance 
quantities. In this method, uncertainties of order a l s remain. Usually, these days, 
only one- loop calculations are available, so / = 2. Then it is helpful to apply the 
Brodsky-Lepage-Mackenzie prescription || to sum up higher-order terms related to 
renormalization parts. The other method is fully non-perturbative [j/J. This sounds 
as if it is exact, but there are uncertainties in the Ks of order a. For the matrix ele- 
ment itself, this is just another error of order a 2 , of which there are many. The first 
method has been used by the QCDSF collaboration ||, who have rather comprehen- 
sive results for the proton. The second method has been used by a Zeuthen-Roma II 
collaboration ||, for the first moment of the pion structure function. 

The quenched results from QCDSF do not agree especially well with phe- 
nomenology. It is, of course, tempting (and reasonable) to blame the quenched 
approximation. Till now, one could also blame the phenomeno logical result, which 
gets a significant contribution from the high-x region, where there are no experimen- 
tal data. As it turns out, neither is the main culprit. Earlier this year the LHPC 
and Sesam collaborations finished a comparison of quenched and unquenched (well, 
n/ = 2) calculations of many proton moments JTTJ. To save computer time, these 
calculations are done with artificially large light quark mass (for the reasons ex- 
plained above). The dependence on the light quark mass for a typical moment is 
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shown in Fig. |l](b). Ref. [HI finds hardly any difference between unquenched and 
quenched calculations for 0.7m s < m q < 1.6m s . It makes a huge difference, how- 
ever, whether one follows one's nose and extrapolates linearly, or whether one follows 
chiral perturbation theory. The latter, of course, is correct. It has a pronounced 
curvature for small quarks masses, of the form lnm^ (and m\ oc m q ). With 
chiral perturbation theory, the extrapolated result agrees with phenomenology. 

3 Flavor Physics 

The central question in flavor physics is whether the standard CKM mechanism 
explains all flavor and CP violation (in the quark sector). One angle on this question 
is over-constraint of the CKM matrix. Because the CKM matrix has only four free 
parameters, the magnitudes of the CKM matrix elements dictate the CP violating 
phase. Many of the magnitudes may be obtained from semi-leptonic decays, such as 
K — > nlv for the Cabibbo angle, B — > D*lv for \Vdb\-, and B — > ixlv for CKM 
elements on the third row (involving the top quark) enter through neutral meson 
mixing, in the neutral K, B, and B s systems. Here we will focus on B physics, with 
a few remarks on K°-K° mixing in Sec. |3.4j . 

For B physics, one must confront heavy quark discretization effects. Com- 
pared to the lattice spacing, the b quark mass is large, m^a > 1. The Symanzik 
effective field theory, at least as usually applied, breaks down. Lattice gauge theory 
does not break down, however, and the Isgur-Wise heavy-quark symmetries emerge, 



in the usual way, for all rriQa [|TTJ. Thus, as long as mg 3> Aq C d, lattice gauge theory 



can be described by heavy-quark effective theory (HQET). One can write |L2| 

Aat = ^C^(mQ,m Q a- it i)O n (ji), (10) 

n 

where = means "has the same matrix elements as." In the same way 

£ QCD = EC nt K;/i)W. (ii) 

n 

The difference is in the short- distance coefficients C n . On the lattice there are two 
short distances, a and mg 1 , so the C^ at depend on the ratio a/nig 1 = rriga. On the 
other hand, the operators O n in Eqs. (|T0|) and ([11]) are essentially the same. 

One can therefore systematically improve lattice calculations of 6-flavored 
hadrons, by matching lattice gauge theory and continuum QCD such that 

5C n = C nt - <t at ~ 0, (12) 
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for the first several operators. HQET is, here, merely an analysis tool; details of 
how HQET is defined and renormalized drop out of the difference. These ideas are 
put to direct use in the Fermilab method fll|, |P2]|, which is based on Wilson fermions 
and, thus, also possesses a smooth continuum limit. Similar ideas are put to use in 
lattice NRQCD IT3], which discretizes the continuum heavy-quark Lagrangian. 



3.1 B -> D*lv, F B ^D'{1), and \V cb \ 

To determine \Va\ from the semi-leptonic decay B — > D*lu, one measures the dif- 
ferential decay rate in w, which is the velocity transfer from the B to the D* . Then, 
one extrapolates to zero recoil, w = 1. Thus, one can summarize the experiment 
by saying it measures \V c b\J-'B-*D*(l), where J 7 b->d*(w) is a certain combination of 
form factors. At zero recoil all form factors but h Al are suppressed, so 

^b-2>.(1) = Ml) = (D*{v)\A"\B(v)). (13) 

It should be "straightforward" to calculate this matrix element in lattice QCD. But 
a brute force calculation of (D*\A tl \B) would not be interesting: similar matrix 
elements like (0|^l M |-B) and (see below) (ix\V^\B) have 15-20% errors. 

At zero recoil heavy-quark symmetry constrains ^^(1) to take the form 

h Al (l) = T] A [l + Si/ m 2 +<5i/ m a] , (14) 

where t\a is a short- distance coefficient, and the 5\/ m n are (principally) long-distance 
matrix elements in HQET at order l/m n . HQET does not provide the tools to 
calculate them, but with the insight from matching lattice gauge theory to HQET, 



we have recently figured out how to do so fL4| . Furthermore, since we incorporate 
heavy-quark symmetry from the outset, and all our uncertainties scale as h Al — 1. 
From HQET, the structure of the 1/mg corrections is 

fi iv i 2iA ip Mtf 

1/m2 (2m c ) 2 + (2m c )(2m b ) (2m 6 ) 2 1 ' 

/3) ,(3)^ ,/3) A /3) 

1/m3 " (2m c )3 + (2m c )(2m 6 ) (2m 6 )3 1 > 

where S = l/(2m c ) + l/(2m&) and A = l/(2m c ) — l/(2m b ). In lattice gauge theory, 
we seek objects whose heavy-quark expansions contain the is. From work on the 
B — ► D form factor [|15|| , we know certain ratios have small enough uncertainties. 



Moreover, one can show via HQET that 

(D\c<y 4 b\B)(B\bj 4 c\D) 
( J D|c7 4 c|D)( J B|67 4 6| J B) " 



{# [l - £ P A 2 - 4 3) A 2 S] } 2 , (17) 



8 



(D*\&f i b\B*) (B*\tr/ c\D*) 
(D*\cj 4 c\D*)(B*\b-f 4 b\B*) 
(D*\c^j 5 b\B)(B*\h^ 5 c\D) 



{Vv 



,A 2 - 



>A 2 



«(3) 



*(3) 



A 2 £ 



A 2 £ 



(18) 
(19) 



{D* | ayj^cl D) {B* | bj^ 5 b\B) 

and one- loop expansions of ify 1 and fj^ are available []16| |. By calculating the ratios 
for many combinations of the heavy-quark masses, we can fit to the HQET descrip- 
tion on the right hand side to obtain all three is in 5i/ m 2, and three of four £^s 
in 5i/ m 3. We can then reconstitute ^^(1) with Eq. ([HP , finding 



n en q+ - 024 -i- n ni«+ - 003 +°- 000 + - 006 
u.y±o_ 0i017 =n u.uiu_ 0i014 _ ^ 016 _ 0>014 , 



(20) 



where the uncertainties stem from statistics and fitting, HQET matching, lattice 
spacing dependence, the chiral extrapolation, and the effect of the quenched ap- 
proximation. Instead of adding these errors in quadrature, we prefer to take note 
of a bound, J-'b^d*^) < 1, and posit a Poisson distribution P(x) = Nx 7 e~ 7x , 
x = [1 - j^ D ,(l)]/0.087 > 0, for global fits of the CKM matrix [0. 



3.2 B -> irlv, f+(E), and \V X 



ub 



To determine \V u b\ from the semi-leptonic decay B — > it is not just a rehash of 
the previous section. The experimental rate is smaller, by a factor \V u b/V c b\ 2 , and 
heavy-quark symmetry is not as constraining. Experiments should measure IT 



dEi 



dE. 



d 2 Y 



n dE t dE T 



oc \V, 



ub 



dE nP 3 \f + {E n ) 



(21) 



where E n = v ■ p n is pion energy in the B rest frame, p 2 = E^ — m 2 , and E\ is 
charged lepton energy. To determine \V u b\ one needs a reliable calculation of the 
form factor f + (E), which parametrizes the matrix element of the b —>■ u vector 
current. Any cut on the lepton variable is equally good 



Recently there have been several calculations of these form factors, using 



several different methods [19, EG, 21, 22, E3y. Two of these works 19, 20 calculate 



the matrix element with rriQ around the charm mass, fit to a model for the E^ 
dependence, and extrapolate the model parameters with heavy-quark scaling. The 
others appeal more directly to heavy-quark ideas, as discussed above. El-Khadra et 
al. [[n[ use the Fermilab method, and the other two |22|, |23| use lattice NRQCD. 



Refs. |21], |22], |23| do not use a model for the E n dependence; instead a cut on E n is 
used to control discretization effects. 

An obvious challenge in these calculations arises from discretization errors 
of the final-state pion, which grow as p n a. This makes it hard to get E n = 2.6 GeV. 
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A possibility to circumvent this difficulty is to give B meson momentum [24]. For 
example, if one chooses — p B = p n = 800 MeV in the lattice frame of reference, 
one can access the whole kinematic range. A less obvious, but also important, 
challenge is the chiral extrapolation. It is not well understood and contributes 



the largest systematic error in the calculation with the smallest quark masses pT [. 
The uncertainties on f+(E) are still 15-20% in the quenched approximation. But 
there are no real technical roadblocks (for details, see Ref. so the errors will 

be reduced while BaBar and Belle accumulate data. In the short term, it will be 
interesting and important to compare similar lattice calculations for semi-leptonic 
D decays to experimental results from CLEO-c. 



3.3 Neutral B Mixing and \V t 



td\ 



In the Standard Model, neutral meson mixing is sensitive to V td and V ts . A signifi- 
cant recent development is the realization that the theoretical uncertainty in B°-B° 
mixing has been underestimated. The culprit has been the chiral extrapolation, 
which we have seen to be important in moments of structure functions. 

In the Standard Model, the oscillation frequency for B®-B® mixing is 

Am d oc \V td \ 2 M d (22) 

where M. q = (S°|[6(l — 7^)759] [&(1 — 7^)75*?] I Phenomenologists usually write 
A4 q = ^m B f B B Bq but lattice calculations give matrix elements M. q directly, and 
f Bq from (0|&7 At 7 5 g|i?g). Nevertheless, it turns out to be useful to look separately at 
fs„ and Bb „■ Current world averages (from lattice QCD) are /r = 198 ± 30 MeV 



and B Bq = 1-30 ± 0.12 [25 . So the error on \Vtd\ from Am^ alone is ~ 15%. 

For some time, the conventional wisdom has said that most of the theo- 
retical uncertainty cancels if one takes the ratio Am s / Am d . (It is anticipated that 
Am, will be measured at Run 2 of the Tevatron E&i.) The ratio is 



Am s 
Am d 



is 



Vt 



id 



m B a ^ c2 _ fB s B B 

fsBs 



e, e = itjf^. (23) 

m B d P b b l 



CKM unitarity says |Vt s | = \V c b\ to good approximation, and \V c b\ is known to 2- 
4%. Many authors believe the uncertainty in £ to be less than 5%. Cancellations 
do occur in the statistical error, and in systematics at short distances (a and m^ 1 ) 
and — arguably — at medium distances (Aq CD ). But they explicitly do not cancel at 
long distances between m~ l and m^ 1 from light quarks in the B s and B mesons. 

Now, the quenched approximation does not work well at these long dis- 
tances, and unquenched calculations are prohibitive at m q ~ m d . Thus, £ isolates 



10 



the contributions that are hardest to capture, and tries to get at them by extrap- 
olating in m q . After studying the differences in chiral logarithms in real QCD and 
the quenched approximation, Booth |27] and Sharpe and Zhang ]2S| sounded notes 
of caution. Their analyses showed that chiral logarithms should induce curvature 
as a function of light quark mass, which quenching would obscure. This curvature 
has recently been observed in unquenched calculations (well, n/ = 2 again), and 
identified as a serious source of uncertainty |2j|. 

It is not too difficult to grasp the problem. For convenience, let £ = £/£_b, 
where = fsj 1 fB d , £s = BbJBb^ Chiral perturbation theory says that 



tf(r)-l =m ss (l-r) 



1 r a- i vy B B*Tr if \ 



(24) 



where m? ss = 2m 2 K — m 2 , r = m q /m s measures the light quark mass in units of the 
strange mass, /„- is the pion decay constant, and gsB*-K is the B-B*-n coupling. The 
function Z(r) contains chiral logarithms: 

(1 - r)Z(r) = |(1 + r) In [(1 + r)/2] + i(2 + r) In [(2 + r)/3] - f r ln(r) 

B s mixing: B s <-> B*K B s <-> B*T] ( 25 ) 

5 d mixing: B «-> S*K B ^ B*r] B «-> B*tt 

and each term arises from the virtual corrections given beneath it. All other con- 
tributions are described well enough by linear behavior in r and are lumped into 
the constant f 2 . The ratio £ B is described by an expression similar to Eq. ([24]), 
except that the chiral log is multiplied by 1 — Qbb*-*- Unquenched lattice calcula- 
tions |29|, ^ are not yet good enough to determine directly the coefficients of the 
chiral logs. Sinead Ryan and I have suggested taking them from phenomenology 
instead JJT]. We invoke heavy-quark symmetry, which says the B-B*-n coupling 
should be roughly the same as the D-D*-tt coupling. Then the recent measurement 
of the D* width yields g 2 DD ^ = 0.35 @; we take g 2 BB . n = 0.35 ± 20%. We obtain 
the constant f'2 from the slope of £/(V) around 0.5 < r < 1.0, where quenched and 
unquenched calculations are in good agreement. We also analyze £ B in the same 
manner. Finally, we find 

£ = 1.32 ± 0.10 (chiral log extrapolation) (26) 

from the chiral log fit. Previously, one had tried linear fits, which would have given 

£ = 1.15 ± 0.05 (conventional linear extrapolation) (27) 

for the same input. The difference is illustrated in Fig. 0(a). One sees that the 
uncertainty of 5% [as in Eq. (27)] certainly was underestimated, and also that the 
central value is probably quite different from the conventional 1.15. 
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3.4 Kaon Physics 



Two theoretical developments of the past few years have opened the way for a wider 
range of kaon calculations. One is a method for exploiting finite-volume effects to 
calculate phase shifts in K — > -nit (and other quasi-elastic processes) . The other 
is the formulation of lattice fermions with good chiral symmetry |34], |3f 



36 . With 



these new tools, we may be able to obtain, at last, quantitative results for such long- 
standing problems as the AJ = ~ rule and the matrix elements needed to compute 
s'/e in the Standard Model. 

Ref. |33[] (and earlier papers by Liischer) base a formalism for calculating 
final-state phases on three insights. The first is that energy levels in finite volumes 
are discrete. The second is that phase shifts arise at hadronic distances of order 
1 fm, remote from the box size L > 1 fm. Lastly, there is a kinematical, albeit 
very non-trivial, set of L-dependent relationships between the phase shifts and the 
energy levels. As a result, one can extract the (elastic) phase shifts from the L 
dependence of the discrete energy spectrum. Unfortunately, practical application to 
QCD requires box sizes L = 2-6 fm, so it will not be feasible soon. 

In many kaonic matrix elements, chiral symmetry is important, for example 
for maintaining a simple relation between lattice operators and their counterparts in 
continuum QCD. In 1982 Ginsparg and Wilson derived a sufficient condition for chi- 
ral symmetry on the lattice |3~7[ . For a long time the Ginsparg- Wilson relation defied 



solution, buy now there are at least two realizations, the "fixed-point action" p5 



and "overlap fermions" ^6[. Another method, "domain- wall fermions" f34|, is re- 
lated to the former and has exponentially small violations of the Ginsparg- Wilson 
relation. Two lattice collaborations have large-scale calculations of Bx (defined 
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Figure 2: Neutral meson mixing, (a) Comparison of chiral extrapolations of £ jffil/ . 
(b) B K vs. a 2 , with results from JLQCD gQj, CP-PACS $ffjjj, and RBC 
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analogously to Bb) with domain- wall quarks. In Fig. |2], published results from 



CP-PACS [[38] and RBC |59] are compared to classic work of JLQCD (with Kogut- 
Susskind quarks) |I(J. The lattice spacing dependence seems gentler for domain- wall 
fermions. From looking at the plot, a rough estimate of an average would be 



B K (NDR, 2 GeV) = 0.58 ± 0.06 



(2? 



which encompasses also Ref. BT 



This uncertainty could easily be reduced, by 
using, say, 3-5 lattice spacings with domain-wall (or overlap) fermions, and taking 
the continuum limit. One should also keep in mind that these calculations have been 
done in the quenched approximation, and with degenerate quarks of mass m s /2. 



4 Conclusions and Prospects 

Although the foundation of lattice QCD is sound, some difficulties arise when turn- 
ing the idealized theory into a computation tool. Errors are introduced at short and 
long distances. They are controlled by effective field theories, however, providing 
reliable methods to obtain physical predictions. A wide variety of calculations in 
the quenched approximation have allowed us to learn how control short-distance 
effects: both for light quarks and for heavy quarks. Now that several (partially) 
unquenched calculations are available, other issues are becoming clearer, particu- 
larly the chiral extrapolation of light quark masses. One can be optimistic that 
unquenched calculations — with solid, transparent analyses of all uncertainties — will 
become available to help interpret experiments with high-energy collisions. 

One upcoming program is especially noteworthy vis a vis lattice QCD. In 
the next few years, CLEO-c [^] will measure leptonic and semi-leptonic decays of 
D and D s mesons to a few per cent. Lattice QCD has a chance to predict their 
results, perhaps with comparable accuracy. An especially interesting combination 
is f^ Klu (E)/f Ds (and the Cabibbo-suppressed cousin f^ nlu (E)/f D ). The CKM 
matrix drops out from the measurements, and non-Standard physics is unlikely. 
Thus, one has direct tests of non-perturbative QCD. The ratio foj fo is also in- 
teresting, because it tests the chiral extrapolation of £ in B-B mixing. Successful 
comparisons of lattice calculations and CLEO-c will give confidence in other appli- 
cations of lattice QCD, such as B physics and moments of the parton densities. 



Acknowledgments 

Fermilab is operated by Universities Research Association Inc., under contract with 
the U.S. Department of Energy. 



13 



References 



1. K. G. Wilson, Phys. Rev. D 10, 2445 (1974). 

2. A. S. Kronfeld, Uses of effective field theory in lattice QCD, in: At the Frontiers 
of Particle Physics: Handbook of QCD, (ed. M. Shifman) 4, (World Scientific, 
Singapore, 2002) [arXiv:hep-lat/0205021]. 

3. W. T. Giele, S. A. Keller, and D. A. Kosower, arXiv:hep-ph/0 104052. 

4. K. Symanzik, Cutoff dependence in lattice 4 in four dimensions theory, in: 
Recent Developments in Gauge Theories, (ed. by G. 't Hooft et al., Cargese, 
August 1979) 59, 313 (Plenum, New York, 1980). 

5. K. Symanzik, Some topics in quantum field theory, in: Mathematical Problems 
in Theoretical Physics, (ed. by R. Schrader et al, Hamburg, August 1981) 6, 47 
(Springer, New York, 1982); Nucl. Phys. B 226, 187, 205 (1983). 

6. S. J. Brodsky, G. P. Lepage, and P. B. Mackenzie, Phys. Rev. D 28, 228 (1983); 
G. P. Lepage and P. B. Mackenzie, Phys. Rev. D 48, 2250 (1993) [arXiv:hep- 
lat/9209022]. 

7. For references and a review, see S. Sint, Nucl. Phys. B Proc. Suppl. 94, 79 (2001) 
[arXiv:hep-lat/0011081], or Ref. 0. 

8. M. Gockeler et al, Phys. Rev. D 53, 2317 (1996) [arXiv:hep-lat/9508004]; Phys. 
Lett. B 414, 340 (1997) [arXiv:hep-ph/9708270]. 

9. M. Guagnelli, K. Jansen, and R. Petronzio, Phys. Lett. B 493, 77 (2000) 
[arXiv:hep-lat/0009006]. 

10. D. Dolgov et al. [LHPC and SESAM Collaborations], Phys. Rev. D 66, 034506 
(2002) [arXiv:hep-lat/0201021]. 

11. A. X. El-Khadra, A. S. Kronfeld, and P. B. Mackenzie, Phys. Rev. D 55, 3933 
(1997) [arXiv:hep-lat/9604004]. 

12. A. S. Kronfeld, Phys. Rev. D 62, 014505 (2000) [arXiv:hep-lat/0002008]. 

13. G. P. Lepage and B. A. Thacker, Nucl. Phys. B Proc. Suppl. 4, 199 (1987); 
B. A. Thacker and G. P. Lepage, Phys. Rev. D 43, 196 (1991); G. P. Lepage, 
L. Magnea, C. Nakhleh, U. Magnea, and K. Hornbostel, ibid. 46, 4052 (1992) 
[arXiv:hep-lat/9205007]. 



14 



14. S. Hashimoto, A. S. Kronfeld, P. B. Mackenzie, S. M. Ryan, and J. N. Simone, 
Phys. Rev. D 66, 014503 (2002) [arXiv:hep-ph/01 10253]. 

15. S. Hashimoto, A. X. El-Khadra, A. S. Kronfeld, P. B. Mackenzie, S. M. Ryan, 
and J. N. Simone, Phys. Rev. D 61, 014502 (2000) [arXiv:hep-ph/9906376]. 

16. J. Harada, S. Hashimoto, A. S. Kronfeld, and T. Onogi, Phys. Rev. D 65, 094514 
(2002) [arXiv:hep-lat/0112045]. 

17. A. S. Kronfeld, P. B. Mackenzie, J. N. Simone, S. Hashimoto, and S. M. Ryan, 
in Flavor Physics and CP Violation, edited by R. G. C. Oldemann, arXiv:hep- 
ph/0207122. 

18. I thank Andreas Warburton for discussions on this point. 

19. K. C. Bowler et al. [UKQCD Collaboration], Phys. Lett. B 486, 111 (2000) 
[arXiv:hep-lat/9911011]. 

20. A. Abada et al, Nucl. Phys. B 619, 565 (2001) [arXiv:hep-lat/0011065]. 

21. A. X. El-Khadra, A. S. Kronfeld, P. B. Mackenzie, S. M. Ryan, and J. N. Simone, 
Phys. Rev. D 64, 014502 (2001) [arXiv:hep-ph/0101023]. 

22. S. Aoki et al. [JLQCD Collaboration], Phys. Rev. D 64, 114505 (2001) 
[arXiv:hep-lat/0 106024]. 

23. J. Shigemitsu, S. Collins, C. T. H. Davies, J. Hein, R. R. Horgan, and G. P. 
Lepage, arXiv:hep-lat/0207011. 

24. K. Foley, poster at the XX th International Symposium on Lattice Field Theory, 
(June 2002, Cambridge, Mass.). 

25. S. Ryan, Nucl. Phys. B Proc. Suppl. 106 (2002) 86 [arXiv:hep-lat/0111010]. 

26. K. Anikeev et al, B physics at the Tevatron: Run II and beyond, arXiv:hep- 
ph/0201071 [FERMILAB-PUB-0 1-197]. 

27. M. J. Booth, Phys. Rev. D 51 (1995) 2338 [arXiv:hep-ph/9411433]. 

28. S. R. Sharpe and Y. Zhang, Phys. Rev. D 53 (1996) 5125 [arXiv:hep- 
lat/9510037]. 



15 



29. N. Yamada et al. [JLQCD Collaboration], Nucl. Phys. B Proc. Suppl. 106 (2002) 
397 [arXiv:hep-lat/01 10087]; talk at the CERN workshop on the CKM Unitarity 
Triangle, (February 2002, Geneva), [http : / /ckm- workshop . web . cern . ch/| . 

30. C. Bernard et al. [MILC Collaboration], arXiv:hep-lat/0206016. 

31. A. S. Kronfeld and S. M. Ryan, Phys. Lett. B 543, 59 (2002) [arXiv:hep- 
ph/0206058]; arXiv:hep-lat/0209083. 

32. A. Anastassov et al. [CLEO Collaboration], Phys. Rev. D 65 (2002) 032003 
[arXiv:hep-ex/0 108043]. 

33. L. Lellouch and M. Liischer, Commun. Math. Phys. 219, 31 (2001) [arXiv:hep- 
lat/0003023]. 

34. D. B. Kaplan, Phys. Lett. B 288, 342 (1992) [arXiv:hep-lat/9206013]; Y. Shamir, 
Nucl. Phys. B 406, 90 (1993) [arXiv:hep-lat/9303005]. 

35. U. J. Wiese, Phys. Lett. B 315, 417 (1993) [arXiv:hep-lat/9306003]; P. Hasen- 
fratz, V. Laliena, and F. Niedermayer, ibid. 427, 125 (1998) [arXiv:hep- 
lat/9801021]. 

36. H. Neuberger, Phys. Lett. B 417, 141 (1998) [arXiv:hep-lat/9707022]; 427, 353 
(1998) [arXiv:hep-lat/9801031]. 

37. P. H. Ginsparg and K. G. Wilson, Phys. Rev. D 25, 2649 (1982). 

38. A. Ali Khan et al. [CP-PACS Collaboration], Phys. Rev. D 64, 114506 (2001) 
[arXiv:hep-lat/0105020]. 

39. T. Blum et al. [RBC Collaboration], arXiv:hep-lat/01 10075. 

40. S. Aoki et al. [JLQCD Collaboration], Phys. Rev. Lett. 80 (1998) 5271 
[arXiv: hep-lat /9 71 00 73] . 

41. G. Kilcup, R. Gupta, and S. R. Sharpe, Phys. Rev. D 57 (1998) 1654 [arXiv:hep- 
lat/9707006]. 

42. R. A. Briere et al, CLEO-c and CESR-c: A New Frontier of Weak and Strong 
Interactions, CLNS-01-1742. 



16 



